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Recently  considerable  efforts  have  been  made  in  understanding  the  dynamics  of 
non-equilibrium  interface  growth  in  the  context  of  a  variety  of  fractal  cluster¬ 
ing  models.  Many  recent  investigations  have  concentrated  on  the  dynamic 
scaling  properties  of  interfaces  obtained  in  experiments  and  in  various  cluster 
growth  models.  Particular  attention  has  been  devoted  to  the  scaling  properties 
of  the  rms  interface  width 

w(l,t)  s  ({k{x,t)  -  ~  (I) 

Here  h{x,i)  is  the  surface  height  at  time  f,  the  angular  brackets  denote  the 
average  over  x  belonging  to  an  interval  of  size  £;  also,  f{u)'^u^  for  ix  ■C  1  and 
f{u)  — »  Const  for  u  ^  1. 

It  has  been  widely  believed  that  many  such  problems  belong  to  the  same 
universality  class  as  the  the  Kardar-Paxisi- Zhang  (KPZ)  equation,® 

?M^=V^h+^{Vhf+^(x,t).  (2) 

Here  T]{x,t)  is  a  random  noise  term.  One  such  model  is  ballistic  deposition,®  for 
which  the  surface  width  exponents  a  and  13  satisfy  the  general  scaling  relation^’® 

a  -t-  ^  =  2  (3a) 

and  can  be  calculated  exactly  in  the  case  of  normally  distributed  uncorrelated 
noise  r}{x,t)  in  one  dimension  (d  =  1  -|- 1): 

«=5.  (36) 

There  have  recently  appeared  several  experiments  on  surface  growth  which  yield 
exponents  quite  different  from  those  of  (3b).  For  example,  in  some  experiments 
"  "  '.iiiniscible  displacement  of  viscous  fluids  in  porous  media  it  was  found  that 
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a  =  0.73±  0.003  [Ref.9],  a  ~  0.81,  /3  =  0.65  [Ref.  10]  and,  in  recent  experiments 
on  the  growth  of  bacteria  colonies  a  =  0.78  ±  0.06  [Ref.llj. 


II.  Roughening  with  Power  Law  Distributed  Noise 


Recently  Zhang^^  suggested  that  the  anomalous  roughening  found  in  experi¬ 
ments  can  be  explained  by  an  uncorrelated  “noise”  T)(x,t)  obeying  a  power-law 
distribution,  p(rf)  ~  where  t)  >  1.  This  anomalous  noise  can  be  simu¬ 

lated  in  a  deposition  context  by  depositing  rods  of  size  i  sampled  from  a  power 
law  probability  distribution, 

(4) 


The  growth  rule  is  similar  to  the  conventional  ballistic  deposition  rule,  i.e.,  a 
deposited  rod  is  attached  to  the  highest  neeurest  neighbor  surface  site.  The  site 
at  which  deposition  next  occurs  can  be  chosen  either  deterministically^^ or 
randomly. 

The  Zhang  idea^^  has  received  recent  support  from  both  mean  field  theory^^ 
and  numerical  simulations,^ both  of  which  suggest  that  the  exponents  a 
and  13  are  anomalously  large,  and  depend  continuously  on  the  parameter  p  (at 
least  for  p  <  pc  ^  5).^^“^® 

Figure  1  shows  the  comparison  of  om  results^"*  for  different  values  of  fi  with 
the  Zhang-Krug  prediction, 


a  = 


fi  +  l' 


2/i-  1’ 


[fi  >  2]. 


(5) 


It  can  be  seen  that  for  /x  >  Pc  =  4.5  ±  0.5  both  a  and  ^  are  almost  independent 
of  /X  Euid  are  very  close  to  the  classical  values  a  =  1/2,  =  1/3.  For  fx  <  pc 

both  exponents  deviate  from  their  “classical”  values  and  approach  the  limiting 
vadues  a  =  1,  /?  =  1  predicted  by  the  Zhang-Krug  relation  for  /x  =  2.  As  a 
final  consistency  check,  we  found  excellent  agreement  with  the  scaling  relation 
Q  -f  q/^  =  2  in  the  entire  range  of  studied  values  of  p  >  2.  For  /x  <  2,  see  our 
discussion  below. 

In  a  recent  work^®  a  closed-form  expression  for  the  probability  distribution 
for  the  fluctuations  8h{x^t)  =  h(x,t)  -  (h(t)}  in  surface  height  h{x,t)  was 
suggested  based  on  a  formal  analogy  between  anomalous  roughening  and  the 
statistics  of  a  Levy  walk. 

By  a  Levy  walk  we  mean  that  each  unit  of  time  a  random  walker  steps  a 
unit  length.  The  walker  moves  £  successive  steps  in  the  same  direction  before 
randomly  changing  direction,  and  £  is  taken  from  a  Levy  distribution^^  p{£)  ~ 
probability  density  P{x,n)  that  the  walker  is  at  position  x  after  n 
steps  has  a  tail  distribution  (x  >■  n*/**)  of  the  form^® 


P{x,n)  ~  n/x'*'*’*  = 


1 


-M-i 


(6) 
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Fig.  1.  Comparison  of  numerical  results  for  exponents  a  (•)  and  (o)  with  theoretical 
predictions  given  for  p  >  2  by  Eqs.  (5)  (solid  line),  and  for  p  <  2  by  (10)  (solid  line). 
After  Ref.  16. 


By  mapping  the  surface  growth  to  Le-^  walk  we  obtain^®  that  for  t  tx  ~ 
L*  the  distribution  of  h(i,  t)  has  the  asymptotic  form 

(7) 

Since  iov  t  ty,,  6h  w  L°,  we  find  a  self  consistent  equation  for  a, 
a  =  (3  —  a)//i,  or  a  =  3/(/i  +  1),  the  same  as  Eq.  (5).  Note  that  a  assumes  its 
classical  value  a  =  1/2  at  /i  =  /ic  =  5.  Substituting  a  =  3/(/i  +  1)  into  (7),  we 
obtain  the  predicted  scaling  form 

(8) 

To  test  (8),  we  performed  simulations  of  ballistic  deposition  using  Elq.  (4) 
for  several  values  of  /i  and  for  a  sequence  of  values  of  L}^  Figure  2  supports 
the  predicted  data  collapse  of  Eq.  (8)  for  /x  =  1.5  and  3. 
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(«)  lool«»ww»  log)5hA»r| 

Fig.  2.  Log-log  scaling  plot  of  wP{6h,L,t)  against  |6/i|/w.  Here  w  is  the  first  moment  of 
P(6h,  L,t)  for  (a),  =  1.5  and  for  (b),  ^  =  3.  Symbols  in  (a)  are:  +  (L  =  1028,  t  =  64);  0 

(L  =  2048,  f  =  256);  x  (X  =  2048,  t  =  1024);  A  (X  =  256,  f  >  4096);  □  (X  =  512,  <  >  4096); 
o  (X  =  1024,  t  >  4096).  Symbols  in  (b)  are;  □  (t  =  32,  X  =  2048);  A  (t  =  128,  X  =  2048);  o 
(t  =  512,  X  =  2048);  -»-(<>  16384,  X  =  512);  x  (t  >  16384,  X  =  1024);  0  (t  >  32768,  X  = 
2048).  The  straight  lines  have  slopes  of  (a)  2.5  and  (b)  4,  as  predicted  by  Eq.  (7).  After  Ref. 
16. 


To  obtain  the  eaxly  time  (l<§:t-Cfx)  dependent  probability  we  use  again 
the  time-spax:e  scaling  relation  t  =  of  (1)  and  L  in  (7)  should  be  replaced 
by  tL  =  yielding 

Since  6h  ~  we  obtain  ^  =  {z-\-\)l fiz  =  3/(2/i  —  1),  the  same  as  in  (5).  Figure 
2  shows  also  data  for  the  time-dependent  probability  density  supporting  (9). 

The  above  considerations,  Eqs.  (7)-(9),  are  valid  for  fi  >  2.  For  fi  <  2, 
the  relation  txL  ^  fails  since  <x  is  bounded  from  below  by  L.  Thus  one 
must  repeat  tl'e  arguments  leading  to  Eqs.  (7)-(9)  using  from  which 

follows^® 


(m  <  2].  (10) 

Note  that  (10)  complements  the  Zhang- Krug  prediction  (5)  for  values  of  p, 
below  2.  Numerical  data  supporting  (10)  are  shown  in  Fig.  1.  The  analogy 


looUoo(P(8  h«0)/P(6H)|| 
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Fig.  3.  (a)  Log-log  scaling  plots  of  log[P(6fc  =  Q,L,t)/P{6h,  L,t)]  versus  \6h\/w  for  p  =  3 
for  different  system  sizes  L  and  times  f .  The  straight  line  has  the  slope  2,  as  should  be  for 
Gaussian  distribution.  The  symbols  are  the  same  as  those  used  in  Fig.  2(b).  (b)  Log-log  plot 
of  log[F(6h  =  0,  L,  t)/P{6h,  L,  t)]  as  a  function  of  \6h\lw  for  different  values  of  p  (the  data  for 
each  value  of  p  obtained  for  L  =  1024  and  large  (,  greater  than  f  x ):  □  (p  =  oo);  x  (p  =  6);  + 
(p  =  5);  0  (p  =  4);  o  (p  =  3);  A  (p  =  2).  The  crossover  value  of  Sh/w  at  which  the  behavior 
changes  from  Gaussian  (straight  line  with  slope  2)  to  power  law  increases  gradually  with  the 
value  of  p.  After  Ref.  16. 


to  Levy  walks  predicts  not  only  the  tails  of  the  probabihty  densities  but  also 
their  behavior  in  the  range  of  small  fluctuations.  In  this  range  and  for  p  >  2 
the  distribution  of  Levy  walks  is  known  to  be  Gaussian,^*  predicting  that  for 
Sh  <  w  the  probability  density  of  stirface  heights  is  also  Gaussian.  Indeed, 
plotting  the  data  in  Fig.  3  as  log{log[P(^h,  L)/P{0,  L)]}  versus  log  6h  for  several 
values  of  L  shows  a  clear  range  of  slope  2  supporting  a  Gaussian  form.  The 
crossover  from  Gaussian  to  a  power-law  occurs  at  a  value  of  y  =  Sh/w  which 
increases  as  p  increases  as  expected  from  the  analysis  of  theory  in  Ref.  18. 

As  seen  from  Fig.  2,  both  time  and  size  dependence  have  the  same  scal¬ 
ing  relation.  Indeed,  Eqs.  (8)  and  (9)  and  the  Gaussian  form  found  for  small 
fluctutions  can  be  combined  for  p  >  2  to  a  single  scaling  relation 

(llo) 
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where  w  =  w{L,t)  is  given  by  Eq.  (1),  and 


F{y)  ~  { 


i  exp(-ay^) 

I 


[y  <  yc] 
[y  >  yc]- 


(life) 


The  data  collapse  shown  in  Fig.  2  supports  (11).  Similar  numerical  results, 
supportng  (11)  for  /i  =  3  were  presented  in  Ref.  13. 

From  the  analogy  to  Levy  walks  we  expect  that  in  the  case  of  conventional 
ballistic  deposition  (/i  =  oo)  the  distribution  will  be  Gaussian; 


(12“) 

for  t  -C  tx,  and 

(126) 

for  t  ^  tx  ~  .  In  Fig.  4  we  show  numerical  data  supporting  (12). 

The  Zhang^^  model,  which  is  based  on  the  assumption  that  the  noise  in  the 
system  has  power  law  distributed  amplitudes,  may  be  the  explanation  for  the 
anomalous  surface  roughening  in  experiments  but  the  origin  of  such  a  noise  in 
real  systems  remains  unclear.^® 


III.  Correlated  Noise 


Anomalous  surface  roughening  can  be  also  due  to  long-range  correlated  noise. 
Next  we  review  recent  studies  on  the  effect  of  long-range  correlated  noise  on 
surface  growth  models. 

When  the  noise  itself  is  the  result  of  another  stochastic  process,  then  the 
noise  cannot  fee  treated  as  random — the  noise  is  correlated  in  space  and/or 
time.^”  In  this  case,  the  exponents  depend  on  the  strength  of  the  correlation. 
Medina  et  al®  used  dynamical  renormalization-group  analysis  to  study  the  KPZ 
equation  with  long  range  correlated  noise.  The  noise  they  studied  has  the  cor¬ 
relation 

~  |x  -  |<  -  (ISa) 

where  d  is  the  overall  dimension  of  the  system  (d  —  1  is  the  dimension  of  the 
surface).  If  the  noise  has  no  temporal  correlation,  i.e., 

(y(x,<)y(x',t'))  ~  -<'),  (13fe) 


the  exponents  obey  the  relation  a  +  z  =  2.  Since  then  there  is  only  one  inde¬ 
pendent  scaling  exponent,  it  is  sufficient  to  give  for  d  =  1  -|-  1 


'1/3 

^il  +  2p)f{5-2p) 


0<P<  1 
J<P<1- 


(14a) 
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Fig.  4:  Log-log  plot  of  log[P(6/i  =  as  a  function  of 

\6h\lui  for  y.  =  oo,  which  is  the  case  of  conventional  ballistic  deposition. 
Here  w  scales  according  to  Eq.  (1)  with  a  =  l/2,;3  =  l/3.  The  Gaussian 
behavior  is  found  in  the  entire  range  of  6h:  □  (f  =  64,  L  =  4096);  A 
(<  =  256,  L  =  4096);  o  {t  =  1024,  L  =  4096);  -f  (L  =  512,  t  >  16384); 
X  (L  =  1024,  <  >  16384);  0  {L  =  2048,  t  >  32768);  •  (L  =  4096,  t  > 
32768).  After  Ref.  16. 


The  other  feature  of  the  KPZ  equation  is  that  it  can  be  mapped  to  the 
directed  polymer  (DP)  problem^^.  The  noise  plays  the  role  of  a  time-dependent 
random  potential.  Thus,  the  results  of  Ref.  5  can  also  apply  to  the  DP  problem 
in  a  correlated  potential  field. 

Zhang^^  used  a  replica  method  to  study  the  DP  problem  with  correlated 
noise  t/  given  by  Eq.  (13b).  Due  to  the  analogy  between  the  DP  problem  and 
the  KPZ  equation,  Zhang  predicts  for  d  =  1  -p  1 


f(l+2p)/(3  +  2^) 
[  (1  +  2/>)/(5  -  ip) 


0  <  p  <  1 

H^<1. 


(146) 
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Hentschel  and  Family^^  studied  the  scaling  behavior  for  dissipative  dynam¬ 
ical  systems  and  proposed  a  new  relation; 


3-2p’ 


o<,<5. 


(14c) 


Note  that  the  three  predictions  [Eqs.  (14a),  (14b),  and  (14c)]  differ  for  0<p< 

1/2. 

There  have  been  several  prior  attempts  to  verify  the  analytical  results  with 
correlated  noise.^^  This  work  relies  on  numerical  methods  that  probably  gener¬ 
ate  undesired  correlations  in  the  noise.  Here  we  review  a  recent  work^®  where  we 
generate  algebraically-correlated  noise,^®  integrate  numerically  the  KPZ  equa¬ 
tion,  and  also  simulate  the  DP  growth  in  a  correlated  potential  field.  The 
results  of  Peng  et  al.^®  for  both  KPZ  and  DP  agree  with  each  other,  and  qual- 
itataively  agree  somewhat  better  with  (14c)  than  with  (14a)  or  (14b).  Finally 
we  implement  correlated  noise  into  the  BD  model,  and  were  surprised  to  find 
surface  roughening  exponents  that  differ  from  both  the  KPZ  equation  and  the 
DP  problem. 

To  construct  the  algebraically- correlated  noise,  we  first  generate  a  represen¬ 
tation  of  random  Gaussian  uncorrelated  noise  T7o(x,t),  then  Fourier  transform 
it  to  obtain  r/o(q,w).  We  define 

^(q,^^)  =  |qr^(‘^("^7o(q,‘*^)-  (is) 


The  noise  r;(x,  t)  is  obtained  by  Fourier  transforming  ^^(q,  a>)  back  into  the  space 
and  time  domain.  It  is  straightforward  to  verify  that  7/(x,  t)  obtained  in  this 
way  has  the  correct  correlations  (13a).  We  restrict  ovurselves  to  the  d  =  1  -|-  1 
case  and  the  noise  has  only  spatial  correlation  (0  =  0)  as  in  Eq.  (13b). 

(i)  Consider  first  the  KPZ  equation  with  noise  rj  described  by  (13b).  For  a 
one- dimensional  surface,  the  discrete  form  of  Eq.  (2)  is 


ht+At{i)  =  ht{i)  +  At[ht{i  +  1)  +  ht{i  -  1)  -  2ht{i)] 


+ 


XAt 


ht{i  +  1)  -  ht{i  -  1) 


l2 


+  '/At'qt{i). 


(16) 


Small  At  is  needed  to  obtain  good  convergence,  and  we  choose  the  appropriate 
time  step  by  verifying  that  smaller  time  steps  do  not  cheinge  our  results.  We 
obtain  the  exponent  /?  from  w{i,i)  defined  in  Eq.  (1),  since  w  t^  for  At  C 
t  <<x. 

We  start  with  the  czise  A  =  0  (no  non-lineaxity)  for  which  z  and  /3  czui 
be  found  exactly  from  dimensional  analysis^ a  change  of  scale  x  bx  and 
t  b^t  implies  h  —*  b^h  and 

7/(x,f) 


[from  Eq.  (13b)].  Equation  (2)  is  scEile  invariant  for  the  choice 


zq  =  2  ; 


(17) 
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Fig.  5.  Comparison  of  our  numerical  results  and  theoretical  predictions  of  (14a),  (14b),  and 
(14c)  (solid,  dashed,  and  dotted  lines  respectively).  Typical  error  bars  are  shown  for  each 
of  the  three  models  treated.  The  dot-dashed  line,  Eq.  (17),  is  obtained  by  neglecting  the 
non-linear  term  in  Eq.  (2).  After  Peng  et  al.,  Ref.  25. 


Our  ntimerical  simulation  for  A  =  0  confirms  (17). 

When  A  0  the  exponents  change.  We  find  that  the  exponent  ^  approaches 
the  same  value  for  non-zero  A.  Since  changing  A  should  not  change  the  tmiver- 
sality  class,  we  czu'ry  out  our  simulation  for  that  value  of  A  which  gives  the 
fzistest  convergence  to  the  correct  value  of  then  we  vary  the  parameter  p. 
The  results  are  shown  in  Fig.  5.  The  solid,  deished,  and  dotted  lines  are  the 
predictions  from  three  theories  [Eqs.(14a),  (14b),  and  (14c)],  respectively. 

To  check  our  results,  we  also  study  the  DP  growth.  By  a  simple  transfor¬ 
mation  W(x,<)  =  exp[(A/2)/i(x,t;j,  we  obtain  from  (2) 

8W  X 

^  +  j,(x,t)W,  (18) 

Here  W  is  the  sum  of  Boltzmann  weights  for  all  configurations  of  a  DP  connect¬ 
ing  (0, 0)  and  (x,  t),  and  fj(x,  t)  is  the  potentieil  field.  The  Boltzmann  weight 
for  all  paths  joining  the  points  (0, 0)  and  (x,  t)  is 

W(x,t)  =  J2expl-Ec/kT].  (19) 

C 

Here  Ec  is  the  sum  of  the  potential  field  rf  on  configuration  c,  and  the  svun  is 
over  all  configvnations  joining  the  two  end  points. 

The  typical  transverse  fluctuation  scales  with  the  length  of  the  polymer  t  as 
(x^(t))^/^  ~  t".  At  zero  temperature,  only  the  optimal  path  (configmration  with 
minimmn  energy)  makes  a  contribution.  Since  the  optimal  path  still  dominates 
at  finite  but  low  temperature,  we  choose  T  =  0  to  simplify  om  numerical  task. 
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We  generate  a  representation  of  [obeying  Eq.  (13b)],  and  record  the  end 

point  of  the  optimal  path  x(t).  We  average  over  many  realizations  (typically  10^ 
of  r]{x,t)).  The  exponent  u  is  related  to  the  dynamic  exponent  z  =  a/ 0  =  2  —  a 
of  KPZ  equation  via  v  =  1/z.  Hence  to  compare  with  the  KPZ  results,  we  define 
ydop  =  2i/  —  1  and  show  the  results  in  Fig.  5.  The  agreement  with  our  numerical 
results  for  the  KPZ  equation  provides  8in  excellent  consistency  check  on  our 
numerical  methods. 

(ii)  Next  we  s^udy  the  BD  model  with  algebreiically  spatial  correlated  noise. 
For  uncorrelated  particles  rain  down  vertically  onto  the  substrate  until 

they  reach  one  of  the  growth  sites.  A  growth  site  is  defined  as  the  highest 
site  on  each  column  that  belongs  to  the  nearest  neighbors  of  the  deposition 
surface.  Once  the  particle  reach  the  growth  site  it  stops  and  become  a  part  of 
the  deposit.  Note  that  the  deposition  rule  defined  above  allows  lateral  growth, 
which  is  believed  to  be  described  by  the  non-linear  term  [(V/i)^]  in  Eq.  (2). 

We  introduce  correlated  noise  according  to  OSb).  As  seen  from  Fig.  5,  we 
find  significant  differences  between  exponents  obtained  from  the  BD  model  and 
the  DP  growth  (or  the  KPZ  model). 


rV.  Experiment  and  a  Directed  Percolation  Model 


Next  we  present  experiments  in  which  ink,  coffee  and  other  suspensions  are 
absorbed  by  a  hanging  paper,  forming  a  rough  interface  between  wet  and  dry 
regions.  We  analyze  this  morphology  and  meeisure  its  roughness  exponent  a,  Eq. 
(1).  Based  on  the  experiment  we  propose  a  new  model  for  interface  roughening. 
Both  the  model  and  the  experiment  produce  interfaces  with  an  anomalously 
large  value  of  a.^* 

(a)  Experiment.  The  experiment  was  performed  by  clipping  paper  to  a  ring 
stand,  and  allowing  it  to  dip  into  a  basin  filled  with  suspensions  of  ink  or  coffee. 
The  suspension  was  absorbed  into  the  paper,  forming  a  rough  interface  between 
the  wet  and  the  dry  regions.  We  allow  the  interface  to  rise  until  it  stops  and 
no  change  in  either  height  or  shape  of  the  interface  is  observed.  The  stopping 
can  be  attributed  to  the  evaporation  of  the  fluid  in  the  wet  regions.  After 
drying,  we  digitize  this  rough  interface.  We  then  calculate  the  height-height 
correlation  function^®  c(£,  0)  =  w(£)  on  different  length  scales  £,  averaging  over 
ten  different  interfaces.  Figure  6a  shows  the  data,  which  support  a  scaling  of 
the  form  c(£,  0)  =  w(£)  ~  £"  with  q  =  0.63  ±  0.04. 

(b)  Model.  The  model  we  propose  is  defined  as  follows:  on  a  square  lattice 
of  edge  L  (with  periodic  boundciry  conditions)  we  block  a  fraction  p  of  the  cells 
to  correspond  to  the  inhomogeneous  nature  of  the  paper  towel.  At  t  =  0,  we 
regard  the  “interface”  to  be  the  bold  horizontal  line  shown  in  Fig.  7a.  At  <  =  1 
we  randomly  choose  a  cell  (labeled  X  in  Fig.  7b)  which  is  one  of  the  unblocked 
dry  cells  that  are  nearest  neighbors  to  the  interface.  We  wet  cell  X  and  any 
cells  that  are  below  it  in  the  same  column.  Ihis  process  is  then  iterated.  For 
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log< 


Fig.  6.  Log-log  plots  shol^  ing  the  dependence  on  length  scale  t  of  the  height-height  cor¬ 
relation  function  c(/,  0)  ~  «;(<)  for  (a)  the  experimental  data  (averaging  over  10  different 
experiments),  and  (b)  the  numerical  results  (averaging  over  1000  different  realizations  for 
system  size  L  =  16, 384  and  for  p  =  0.4675,  very  close  to  pc  for  the  infinite  system).  The 
slope  for  the  set  of  experimental  points  indicated  by  solid  circles  (two  decades)  is  0.63  ±0.04, 
while  the  slope  for  the  simulation  point  indicated  by  solid  circles  (three  decades)  is  0.63±0.02. 
After  Barabasi  et  al.  Ref.  28. 


example,  Fig.  7c  shows  that  at  t  =  2  we  choose  cell  Y  a  second  unblocked  cell 
to  wet,  while  Fig.  7d  shows  that  at  t  =  3  we  wet  cell  Z  and  also  cell  Z'  below 
it?^ 

We  find  that  for  p  below  a  critical  threshold^®  Pc  =  Pc{L)  the  interface 
propagates  without  stopping,  while  for  p  above  pc  the  interface  is  pinned.  Figure 
6b  displays  the  scaling  behavior  of  the  model  at  criticality,  and  we  fiind  that 
o  =  0.63  db  0.02,  a  value  identical  to  the  experimentsJ  value  of  Fig.  6a. 
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Fig.  7.  Explanation  of  the  model  for  interface  growth  with  erosion  of  overhangs.  Wet  cells 
are  indicated  by  shaded  cells.  Dry  cells  are  randomly  blocked  with  probability  p  (indicated 
by  0  or  unblocked  with  probability  1  —  p  (indicated  by  1).  The  interface  between  wet  and  dry 
cells  are  shown  by  a  heavy  line,  (a)  <  =  0,  (b)  <  =  1,  (c)  t  =  2,  and  (d)  <  =  3.  After  Barabasi 
et  al.,  Ref.  28. 


Next  we  eirgue  that  the  model  presented  above  is  connected  to  directed 
percolation,^ thereby  providing  a  theoretical  basis  for  the  observed  emd  cal¬ 
culated  values  of  the  anomalous  roughening  exponent  q.  The  propagation  of  the 
interface  will  stop  when  it  reaches  for  the  first  time  a  directed  path  of  blocked 
cells  leading  from  West  to  East — this  path  is  such  that  one  can  walk  on  it  from 
West  to  Eeist  without  turning  to  the  West.  Such  a  ‘directed  path’  is  a  path 
on  the  directed  percolation  cluster  formed  by  the  cells  labeled  0.  We  assume 
that  a  single  transverse  length  characterizes  the  directed  percolation  clusters 
so  that  the  width  w  of  this  interface  scales  as  the  transverse  correlation  length 
of  the  directed  percolation  problem  (^x  is  a  rigorous  upper  boimd).  Thus 
we  assume  w(£)  ~  ^x  and  £  ~  ^K,  where  ^jj  is  the  longitudinal  correlation 

•  •  I 

length  in  the  corresponding  directed  percolation  problem.  Since  ^x  ~  ^|| 
we  identify^ a  =  t'n/i'x  —  0.63. 

To  probe  the  dynamics  of  the  growing  interface  in  the  model,  we  study 
the  height-height  correlation  function  cf£,t).  Our  niimerical  results  support 
an  exponent  ^  =  0.64  ±  0.04.  The  usual  exponent  identity  attributed  to  the 
Galilean  invariance,  which  is  known  to  be  valid  for  the  KPZ  equation,®  is  vi¬ 
olated;  we  find  a  +  2  Cii  1.69,  smaller  than  two.  This  is  a  consequence  of  the 
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strong  anisotropy  of  the  mechanism  which  excludes  the  overhangs:  An  infinites¬ 
imal  tilting  of  the  pinned  interface  will  result  in  removing  blocked  cells,  thus 
allowing  the  interface  to  propagate  further. 

Further  support  for  the  directed  percolation  model  can  be  obtEuned  if  we 
consider  a  finite  system  at  a  fixed  value  po  <  Pc-  If  ^||(Po)  is  larger  than  the 
system  size  L,  the  interface  may  be  stopped  by  the  directed  percolation  path. 
Thus  we  identify  two  regimes:  Regime  I  where  ^||  >  L  and  Regime  II  where 
^11  <  L.  In  Regime  I,  we  observe  only  anomalous  roughening  (a  ~  0.68),  while 
in  Regime  II  we  predict  a  crossover  to  behavior  described  by  the  KPZ  exponent 
(a  =  0.5).  A  similar  crossover^^  is  observed,  both  in  our  calculations  and  even 
in  some  recent  experiments  (Fig.  3  of  Ref.  10). 


Fig.  8.  A  two-dimensional  “directed  surface”  with  erosion  of  overhangs  using  the  rules  de¬ 
scribed  in  Fig.  7.  The  critical  probability  of  blocked  sites  is  pc  —  0.74. 


We  also  studied  several  variants  of  the  above  model.  One  interesting  vari¬ 
ant  arises  if  we  replace  the  blocked-unblocked  percolation  substrate  by  one  used 
in  invasion  percolation.  At  every  cell  of  the  lattice  we  put  a  random  number 
between  0  and  1,  and  advance  the  interface  to  the  nearest  neighbor  with  the 
smallest  random  number.  We  use  the  same  mechanism  as  in  the  previous  model 
to  erode  the  overhangs.  In  this  model  the  interface  never  stops  propagating  im- 
less  we  introduce  a  cutoff  (a  particular  value  of  p  =  pi ;  cells  possessing  a  random 
number  p>  p\  are  blocked  cells).  Thus  the  interface  that  stops  propagating  is 
generated  by  the  same  mechanism  as  in  the  normal  percolation  model,  i.e.  they 
are  in  the  same  universality  class.  However,  one  significant  difference  is  that 
the  normal  percolation  interface  is  generated  basically  by  a  local  growing  rule, 
while  in  the  invasion  percolation  the  growing  rule  is  global.  Numerical  studies 
on  the  moving  interf2w:e — both  of  invasion  and  normal  percolation  substrate — 
give  the  same  roughness  exponent,  a  =  0.68  ±  0.05,  slightly  larger  than  for  the 
pinned  interface. 

Barabasi  et  al  also  studied  the  above  model  in  higher  dimensions,  for 
which  there  are  no  theoretical  predictions  on  the  values  of  the  scaling  exponents 
(see  Fig.  8).  Our  results  suggest  a  =  0.51  ±  0.05  for  2  H-  1  dimensions,  which 
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is  larger  than  the  most  accurate  available  numerical  result^^  for  KPZ  growth 
a  =  0.4  ±  0.01.  Using  similar  arguments  as  for  =  1  +  1,  we  obtain  here  a  two- 
dimensional  “directed  surface”  with  correlation  exponents  u\\  and  which 
obey  a  =  u\\/u±  =  0.51. 
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of  Disorderly  Growth:  Recent  Approaches  to  Understanding  Diffusion  Limited  Aggregation,” 
Physica  A  168,  23-48  (1990).  [invited  talk,  international  conference  on  frontiers  in 

CONDENSED  MATTER  PHYSICS.] 

7.  S.  Havlin,  A.-L.  Barabasi,  S.  V.  Buldyrev,  C.  K.  Peng,  M.  Schwartz,  H.  E.  Stanley,  and  T. 
Vicsek,  “Anomalous  Surface  Roughening:  Experiment  and  Models,”  in  Growth  Patterns  in 
Physical  Sciences  and  Biology  [proc  i99i  NATO  advanced  research  workshop,  Granada, 
Spain,  October  1991],  E.  Louis,  L.  Sander  and  P.  Meakin,  eds.  (Plenum,  New  York,  1992). 

8.  A.  Coniglio,  “Correlations  in  Thermal  and  Geometrical  Systems,”  in  Correlations  and  Con¬ 
nectivity:  Geometric  Aspects  of  Physics,  Chemistry  and  Biology  [NATO  asi  series  vol.  iss],  eds. 
H.  E.  Stanley  and  N.  Ostrowsky  (Kluwer,  Dordrecht,  1990). 

RESEARCH  ARTICLES  PUBLISHED  IN  REFEREED  JOURNALS  DURING  THE  GRANT  PERIOD 

9.  P.  Alstrpm  and  F.  Sciortino,  “Dynamics  of  Bonded  Networks  with  Two  Energy  Scales,”  Phys. 
Rev.  Lett.  65,  2885  (1990). 

10.  P.  Alstrpm,  P.  Trunfio,  and  H.  E.  Stanley,  “Spatio-Temporal  Fluctuations  in  Growth  Phe¬ 
nomena:  Dynamical  Phases  and  1/f  Noise,”  Phys.  Rev.  A  Rapid  Communications  41,  3403 
(1990). 

11.  F.  Caserta,  H.  E.  Stanley,  W.  Eldred,  G.  Daccord,  R.  Hausman,  and  J.  Nittmann,  “Physical 
Mechanisms  Underlying  Neurite  Outgrowth:  A  Quantitative  Analysis  of  Neuronal  Shape,” 
Phys.  Rev.  Lett.  64,  95-98  (1990). 

12.  P.  Devillard  and  H.  E.  Stanley,  “Exact  Enumeration  Approach  to  the  Directed  Polymer  Prob¬ 
lem,”  Phys.  Rev.  A  41,  2942-2951  (1990). 

13.  J.  Lee,  A.  Coniglio,  and  H.  E.  Stanley,  “Fractal-to-Nonfractal  Crossover  for  Viscous  Fingers,” 
Phys.  Rev.  A  Rapid  Communications  41,  4589-4592  (1990). 

14.  J.  Lee,  S.  Havlin,  H.  E.  Stanley,  and  J.  E.  Kiefer,  “Hierarchical  Model  for  the  Multifractality 
of  Diffusion  Limited  Aggregation,”  Phys.  Rev.  A  42,  4832-4837  (1990) 

15.  T.  Nagatani  and  H.  E.  Stanley,  “Double  Crossover  Phenomena  in  Laplacian  Growth:  Effects 
of  Sticking  Probability  and  Finite  Viscosity  Ratio,”  Phys.  Rev.  A  41,  3263-3269  (1990). 

16.  T.  Nagatani  and  H.  E.  Stanley,  “Phase  Transition  and  Crossover  in  Diffusion-Limited  Aggre¬ 
gation  with  Reaction  Times,”  Phys.  Rev.  A  42,  3512-3517  (1990). 

17.  T.  Nagatani  and  H.  E.  Stanley,  “Crossover  and  Thermodynamic  Representation  in  the  Ex¬ 
tended  T]  Model  for  Fractal  Growth,”  Phys.  Rev.  A  42,  4838-4844  (1990). 

18.  C.-K.  Perm,  S.  Prakash,  H.  J.  Herrmann,  and  H.  E.  Stanley,  “Randomness  versus  Deterministic 
Chaos:  Effect  on  Invasion  Percolation  Clusters,”  Phys.  Rev.  A  42,  4537-4542  (1990). 

19.  S.  Schwarzer,  J.  Lee,  A.  Bunde,  S.  Havlin,  H.  E.  Roman,  and  H.  E.  Stanley,  “Minimum  Growth 
Probability  of  Diffusion- Limited  Aggregates,”  Phys.  Rev.  Lett.  65,  603-606  (1990). 
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20.  F.  Sciortino,  A.  Geiger,  and  H.  E.  Stanley,  “Isochoric  Differential  Scattering  Functions  in  Liquid 
Water:  The  Fifth  Neighbor  as  a  Network  Defect,”  Phys.  Rev.  Lett.  65,  3452  (1990). 

21.  F.  Sciortino,  P.  Poole,  H.  E.  Stanley  and  S.  Havlin,  “Lifetime  of  the  Hydrogen  Bond  Network 
and  Gel-Like  Anomalies  in  Supercooled  Water,”  Phys.  Rev.  Lett.  64,  1686-1689  (1990). 

22.  R.  L.  B.  Selinger  and  H.  E.  Stanley,  “Percolation  of  Interacting  Diffusing  Particles,”  Phys. 
Rev.  A  42,  4845-4852  (1990). 

23.  H.  E.  Stanley,  P.  Poole,  A.  Coniglio,  and  N.  Jan,  “What  is  the  Shape  of  a  Polymer  Chain  near 
the  Theta  Point?”  Mol.  Cryst.  &  Liquid  Crst.  18,  91-92  (1990). 

24.  D.  Stassinopoulos,  G.  Hnber,  and  P.  Alstrpm,  “Measuring  the  Onset  of  Spatio-Temporal  In- 
termittency,”  Phys.  Rev.  Lett.  64,  3007-3010  (1990). 

25.  B.  Stosic,  S.  Milosevic,  and  H.  E.  Stanley,  “Density  of  States  Functions  of  the  FuUy  Finite 
Ising  Model  Systems,”  Phys.  Rev.  B  41,  11466  (1990). 

26.  P.  Trunfio  and  P.  Alstrpm  “Exponentially  Small  Growth  Probabilities  in  Diffusion- Limited 
Aggregation,”  Phys.  Rev.  B  41,  896-898  (1990). 

27.  M.  Araujo,  S.  Havlin,  and  H.  E.  Stanley,  “Anomalous  Fluctuations  in  Tracer  Concentration  in 
Stratified  Media  with  Random  Velocity  Fields,”  Phys.  Rev.  A  44,  6913  (1991). 

28.  M.  Araujo,  S.  Havlin,  G.  H.  Weiss,  and  H.  E.  Stanley,  “Diffusion  of  Walkers  with  Persistent 
Velocities,”  Phys.  Rev.  A  43,  5207-5213  (1991). 

29.  M.  Araujo,  S.  Havlin,  H.  Larralde  and  H.  E.  Stanley,  “‘Random- Force  Dominated’  Reaction 
Kinetics:  Particles  Moving  Under  Random  Forces,”  Phys.  Rev.  Lett.  68,  1791  (1992) 

30.  A.-L.  Barabasi,  S.  V.  Buldyrev,  F.  Caserta,  S.  Havlin,  H.  E.  Stanley,  and  T.  Vicsek,  “Anomalous 
Interface  Roughening  in  Porous  Media:  Experiment  and  Model,”  Phys.  Rev.  Lett,  (submitted). 

31.  K.  R.  Bhaskar,  D.  Gong,  B.  S.  Turner,  R.  Bansil,  H.  E.  Stanley,  and  J.  T.  LaMont,  “Viscous 
Fingering,  a  Physico-Chemical  Property  of  Gastric  Mucin,  May  Contribute  to  Gastric  Mucosal 
Acid  Resistance,”  Gastroenterology  100,  A34  (1991). 

32.  K.  R.  Bhaskar,  B.  S.  Turner,  P.  Garik,  J.  D.  Bradley,  R.  Bansil,  H.  E.  Stanley,  and  J.  T. 
LaMont,  “  Viscous  Fingering  of  HCl  Through  Gastric  Mucin,”  Nature  (submitted). 

33.  S.  V.  Buldyrev,  S.  Havlin,  J.  Kertesz,  H.  E.  Stanley,  and  T.  Vicsek,  “Ballistic  Deposition  with 
Power  Law  Noise:  A  Variant  of  the  Zhang  Model,”  Phys.  Rev.  A  Rapid  Commun.  43,  7113 
(1991). 

34.  H.  Larralde,  P.  Trunfio,  S.  Havlin,  H.  E.  Stanley,  and  G.  H.  Weiss,  “Territory  Covered  by  N 
Diffusing  Particles,”  Nature  355,  423  (1992), 

35.  H.  Larralde,  M.  Araujo,  S.  Havhn,  and  H.  E.  Stanley,  “Reaction  Front  for  A  +  B  C  Diffusion 
Reaction  Systems  with  Initially  Separated  Reactants,”  Phys.  Rev.  A  (submitted). 

36.  H.  Larralde,  P.  Trunfio,  S.  Havlin,  H.  E.  Stanley,  and  G.  H.  Weiss,  “Number  of  Distinct  Sites 
Visited  by  N  Random  Walkers,”  Phys.  Rev.  A  45,  xxx  (1992). 

37.  J.  Lee,  S.  Havlin  and  H.  E.  Stanley,  “Analytic  Solution  of  the  Growth  Site  Probability  Dis¬ 
tribution  for  Structural  Models  of  Diffusion  Limited  Aggregation,”  Phys.  Rev.  A  45,  1035 
(1992). 

38.  T.  Nagatani  and  H.  E.  Stanley,  “Viscous  Fingering  near  Percolation  Threshold:  Double- Crossover 
Phenomena,”  Phys.  Rev.  A  43,  2963-2969  (l991). 

39.  T.  Nagatani,  J.  Lee,  and  H.  E.  Stanley,  “Renormalization  Group  for  Viscous  Fingering  with 
Chemical  Dissolution,”  Phys.  Rev.  Lett.  66,  616  (1991). 

40.  T.  Nagatani,  J.  Lee,  and  H.  E.  Stanley,  “Crossover  Effects  in  Chemical  Dissolution  Phenomena: 

A  Renormalization  Group  Study,”  Phys.  Rev.  A  45,  2471  (1992). 

41.  C.  K.  Peng,  S.  Havlin,  M.  Schwartz,  and  H.  E.  Stanley,  “Directed-Polymer  and  Ballistic- 
Deposition  Growth  with  Correlated  Noise,”  Phys.  Rev.  A  44,  2239-2242  (1991). 

42.  C.  K.  Peng,  S.  Buldyrev,  A.  Goldberger,  S.  Havlin,  F.  Sciortino,  M.  Simons  and  H.  E.  Stanley, 
“Long-Range  Correlations  in  Nucleotide  Sequences,”  Nature  355,  168-171  (1992). 

43.  S.  Sastry,  F.  Sciortino,  and  H.  E.  Stanley,  “Collective  Excitations  in  Liquid  Water  at  Low 
Frequency  and  Large  Wave  Vector,”  J.  Chem.  Phys.  95,  7775  (1991). 
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44.  S.  Schwarzer,  J.  Lee,  S.  Havlin,  H.  E.  Stanley,  and  P.  Meakin,  “Distribution  of  Growth  Prob¬ 
abilities  for  Off-Lattice  Diffusion- Limited  Aggregation,”  Phys.  Rev.  A  43,  1134-1137  (1991). 

45.  F.  Sciortino,  R.  Bansil,  H.  E.  Stanley,  and  P.  Alstrom,  “Interference  of  Phase  Separation  and 
Gelation,”  Phys.  Rev.  Lett,  (submitted). 

46.  F.  Sciortino,  A.  Geiger,  and  H.  E.  Stanley,  “Effect  of  Defects  on  Molecular  Mobility  in  Liquid 
Water,”  Nature  354,  218  (1991). 

47.  F.  Sciortino,  A.  Geiger,  and  H.  E.  Stanley,  “Network  Defects  and  Molecular  Mobility  in  Liquid 
Water,”  J.  Chem.  Phys.  96,  3857  (1992). 

48.  C.  Amitrano,  A.  Coniglio,  P.  Meakin  and  M.  Zannetti  “Multiscaling  in  Diffusion  Limited 
Aggregation,”  Phys.  Rev.  B  44,  4974  (1991). 

49.  A.  Coniglio,  F.  di  Liberto,  G.  Monroy  and  F.  Peruggi,  “Cluster  Approach  to  Spin  Glasses  and 
the  Frustrated  Percolation  Problem,”  Phys.  Rev.  B  44,  12,605  (1991). 

50.  A.  Coniglio,  Y.  Oono,  A.  Shinozaki  and  M.  Zannetti,  “On  the  Nature  of  Dynamic  Scaling  in 
Spinodal  Decomposition,”  Europhys.  Lett.  18,  59  (1992). 

51.  S.  C.  Glotzer,  F.  Sciortino,  M.  Gyure,  R.  Bansil,  and  H.  E.  Stanley,  “Microphase  Separation 
in  Polymers:  Interference  between  Thermoreversible  Gelation  and  Spinodal  Decomposition,” 
Nature  (submitted). 

52.  J.  Lee,  A.  Coniglio,  and  H.  E.  Stanley,  “Localization  of  Growth  Sites  in  DLA  Clusters:  Multi- 
fractality  and  Multiscaling,”  Nature  (submitted). 

53.  C.  K.  Peng,  A.  L.  Goldberger,  J.  Hausdorff,  J.  Mietus,  and  H.  E.  Stanley,  “Beat-to-Beat  Heart 
Rate  Fluctuations:  Long-Range  Correlation  and  Their  Breakdown  in  Cardiac  Disease,”  Nature 
(submitted). 

54.  F.  Caserta,  R.  E.  Hausman,  W.  D.  Eldred,  H.  E.  Stanley,  and  C.  Kimmel,  “Effect  of  Viscosity 
on  Neurite  Outgrowth  and  Fractal  Dimension,”  Neurosci.  Letters  136  198-202  (1992). 

55.  F.  Sciortino,  H.  E.  Stanley,  J.  Shao,  G.  Wolf,  and  C.  A.  AngeU,  “Stability  Limits  at  Positive 
and  Negative  Pressures  for  Tetrahedrally- Coordinated  Crystals,”  Phys.  Rev.  Lett,  (semifinal 
draft). 


oiir-annual5.tex:  List  of  awards,  invited  talks,  and  publications 


Page  4 


